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ABSTRACT 

Recently, strong AGN outbursts at the centers of galaxy clusters have been 
found. Using a simple model, we study particle acceleration around a shock 
excited by an outburst and estimate nonthermal emission from the accelerated 
particles. We show that emission from secondary electrons is consistent with the 
radio observations of the minihalo in the Perseus cluster, if there was a strong 
AGN outburst > 10^ yrs ago with an energy of ~ 1.8 x 10®^ erg. The validity 
of our model depends on the frequency of the large outbursts. We also estimate 
gamma-ray emission from the accelerated particles and show that it could be 
detected with GLAST. 

Subject headings: acceleration of particles — radiation mechanisms: non-thermal 
- galaxies: active — galaxies: clusters: general — galaxies: clusters: individual: 
Perseus (A426) 



1. Introduction 



Diffuse synchrotro n radio ernission is often found in the i ntracluster medium (ICM] 



of galaxy clusters (e.g., 


Kim et al. 


1990: 


Giovannini et al. 


1993; 


Giovannini & Feretti 


2000; 


Kempner & Sarazin 


2001 


). These radio sources in clusters are often classified as either 
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peripheral cluster radio relic sources or central cluster radio halo sources. Since most 
of the relics and halos are found in merging clusters, the particles responsible for the 
emission are thought to be accelerated at sho c ks or turbulence excited by cluster mergers 

Takizawa fc N^itolboool: 



(iRoettiger et al 



19991: 



Fuiita fc Sarazin 



Ohno. Takizawa. fc Shibatal I2OO2I : iFujita. Takizawa. fc SarazinI l2003l : iBrunetti et al.l 120041 ) . 



2001 



Brunetti et al 



2001 



However, there are exceptional diffuse radio sources. They are called 'mi nihalos' and are 



1991; 


Burns et al. 


1992; 


Rizza et al. 


2000 


)• 


Gitti et al. 


(2002) 



chrotron emission from radio minihalos is due to a relic population of relativistic electrons 
reaccelerated by MHD turbulence via secon d-order Fermi acceler a tion, and that the energy 
is supplied by cooling flows. Alternatively, iPfrommer &: EnBlinI (120041 ) discussed that the 
electrons responsible for the synchrotron emission from minihalos are o f secondary origin 
and thus are injected during proton-proton collision in t he ICM. However. iGitti et al.l (120021 ) 



did not discuss the generation of the turbulence, and iPfrommer fc EnBlinI (120041 ) did not 
investigate the mechanism of proton acceleration. 

Recent X-ray observations have shown that AGNs at the centers of galaxy clusters 
sometimes exhibit inte nsive outbursts with a mechanical power of ~ 10^^ erg. The examples 



are M S 0735.6-^7421 (IMcNamara et al 



2005 



Gitti et al.l 120071) . Hercules A (iNulsen et al. 



2005ah . and Hydra- A jNulsen et al.l l2005bl : IWise et al.l 120071 ). Such an intensive outburst 



should excite a shock in the ICM. In fact, weak shocks have been found in those clusters. 
In the early stage of the evolution of the shock, the Mach number is expected to be large. 
Therefore, particles would be accelerated at the shock as in the case of a supernova remnant. 

In this letter, we consider particle acceleration at the shock generated by an intensive 
outburst of the AGN at the cluster center. We study nonthermal emission from the accel- 
erated pa rticles. In particular, we fo cus on nonthermal radio emission of secondary origin. 
Recently, iHinton fc Domainkd (120071 ) estimated gamma-ray emission associated with AGN 
outbursts, assuming that the hot cavity behind a shock is entirely filled with high-energy 
protons. This assumption may be too simple, and they did not discuss the acceleration and 
energy spectrum of protons. However, motivated by this study, we also estimate the gamma- 
ray emission using our model. We take the Perseus cluster as a model cluster, because this 
cluster has a well-studied minihalo. 



3 



Models 



We assume that the duration of an AGN outburst is ~ 10'' yr. Since we are interested 
in the evolution of a shock for t > 10^ yr, where t = corresponds to the ignition of the 
outburst, the evolution can be approximated by that of an instant explosion with an energy 
of at t = 0. 

For the sake of simplicity, we assume that the cluster is spherically symmetric and the 
density profile of the ICM has a form of a power-law: 



picuir) = piir/riY 



-3 



10 kpc, and u = 1.43, based on the den sity profile 



We take pi = 5.3 x 10"^^ g cm~", ri 
of the Perseus cluster for 70 < r < 300 kpc (Fig. 8 in IChurazov et al.l 120031 ). We take 
that region because our model is correct for t > 1 x 10'' yr, and the radius of the shock at 
t ~ 10'' yr is Rs ~ 70 kpc for parameters we adopt in ^ 



Using a shell approximation (e.g. 
be written as 



Ostriker fc McKedll988l ). the radius of the shock can 



pirt 



where 



3 (7 + l)2(7-l)(3-cj) 
An 97-3-cj(7+l) 



l/(5-c^) 



(2) 



(3) 



and 7(= 5/3) is the adiabatic index. The Mach number of the shock gradually decreases. The 
shock stops expanding when its Mach number approaches to one. At this point, the cavity 
filled with hot gas inside the shock becomes in pressure equilibrium with the surrounding 
ICM. Since we stop calculation before the radiative cooling of the shock becomes effective 
(> Gyr), we do not need to consider the radiative cooling. 



Following lYamazaki et al.l (120061 ). we assume that particles are accelerated at the shock 
via diffusive shock acceleration (i.e., first-order Fermi acceleration) and that their energy 
spectra are given by 



N{E) oc E-^'e 



-E/E^ 



(4) 



where -E'max is the maximum energy of the protons or electrons. T he index is given by 
X = (rh + 2)/{rh — 1), where rj, is the compression ratio of the shock (IBlandford fc Eichler 



19871 ). We estimate the maximum energies of the protons and electrons using the relations 
of 



^acc mill{^pp) t} ) ^acc mill{^syn) t} ; 



(5) 
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respectively. Here, tacc, tpp, t and tgyn are the acceleration time, the hfetime of high-energy 
protons through pion production, the age of the shock wave, and the synchrotron coohng 
time, respectively. 

Assuming the standard manner for the diffusion coefficient, the acceleration time is 
given by 

_ 20/lC-5'niax 

where c is the speed of light, — e is the el ectron charge, and K(= dRs/dt) is the shock 
velocity (jJokipiil 119871 : lYamazaki et al.ll2004l ). The correction factor h depends on the mean 
free path of particles and the angle between the shock and the magnetic field. Since /i ~ 1 
in the Bohm-limit case, we assume that h = 1 from now on. The downstream magnetic field 
is given by Bd = r^B, where B is the magnetic field strength of the unperturbed ICM. We 
estimate tpp as 



tpp = 5.3 X lO''' yr (nicM/cm 



(7) 



where njcu is the number density of the ICM. Since the shock is in pressure equilibrium in 
~ 10® yr (see §3]) and riicu ^ 0.1 cm"^, the cooling is not effective. Thus, the maximum 
energy of protons is determined by the age of the shock. On the other hand, the synchrotron 
cooling time for electrons is given by 



t 



syn 



1.25 X lOVr 



-1 



B, 



-2 



lOTeVy VlO/iG, 
and is shorter than the age of the shock. Thus, using relations we obtain 

K V f Bd \ f t 



E ~ 



1.6 X 10' 



E ~ 



14 



10^ km s ^ 

Vs 

10^ km s~^ 



10 12 G J VlO^yr 

Bd ^-^/^ 



TeV 



TeV. 



(9) 
(10) 



We assume that the minimum electron and proton energies are their rest masses. For 
given proton and electron spectra, we calculate radiation from them. We consider the syn- 
chrotron, bremsstrahlung, and inverse Compton emissions from primary electrons that are 
directly accelerated at the shock, the 7r°-decay gamma-ray through proton-proton collisions, 
and the synchrotron, bremsstrahlung, and inverse Compton emissions from secondary elec- 
trons created thro ugh the decay of charged pions that are also generate d through proton- 
proton collisions (ISturner et al.l 119971 : iKohri. Yamazaki. &: Bambal 120071 ) . The density of 



target protons for the proton-proton interaction is given by r5picM(-Rs)/(l-4 rrip), where rrip 
is the proton mass. We assume that the spectrum from secondary electrons is stationary 
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if the lifetime of the electrons is smaller than the age of the system. On the other hand, 
if the lifetime is larger than the syst em age, we calcul ate the evolution according to §3 of 



Atoyan fc Aharoniaru (119991 ) (see also lKohri et al.ll2007l ). 



Results 



In our model, the evolution of a shock is determined by picM and Eq (equation [2]). The 
Mach number also depends on the ICM temperature, T. The energy spectrum of particles 
depends on the evolution of the shock and the magnetic field, B. The luminosity of nonther- 
mal emission from the shock depends on the total energy of high-energy (> mpC^) protons 
inside the shock, cEq, where < e < 1. We fix Picm('^), T, and B{r) from observations. On 
the other hand, we regard Eo and e as fitting parameters, because there are no observational 
data for them. 

We assume tha t T = 3.5 keV, whi ch is the temperature of the central region of 
the Perseus cluster (IChurazov et al.ll2003l ). although the temperature before the outburst 
might have been somewhat lower. As far as we know, there are no observations of mag- 
netic fields in the Perseus cluster at r > 70 kpc. On the other hand, the observations 
of Faraday rotation s howed that the typical magnetic fiel d strength in clusters for r < 
500 kpc is 5-10 fiG (jClarke. Kronberg. fc Bohringerl 1200 ll ). Therefore, we take B{r) = 
IfiG (picM[^]/picM[150kpc])^/^ assuming that the magnetic field is adiabatically compressed. 
We note that the spectra of particles (equations [9] and [10]) and synchrotron emission from 
high-energy electrons depend on B. Thus, the results shown below is fairly sensitive to the 
assumption on B. 

In the following, the energy of an AGN explosion is En = 1.8 x 10^^ erg, whic h is three 



times larger than the one observed for MS 0735.6-1-7421 (IMcNamara et al.ll2005l ). We use 
this value to match Rg with the size of the radio minihalo in the Perseus cluster. We take 
the acceleration efficiency of e = 0.05 to match radio observations (see below). The ratio 
of high-energy electrons to high-energy protons is taken to be r^-p = 1/1000 as a fiducial 
value. In the self-similar solution of the shock we adopt, the kinetic and thermal energies 
are respectively constant. Therefore, we assume that the total energy of the high-energy 
protons is also constant. The typical Mach number of the shock in our calculations is ~ 3 
(for 1 X 10'' < t < 4 X 10'' yr). Performing simulations takin g account of the back-reaction 
of accelerated particles on hydrodynamics, iRyu et al.l (120031 ) estimated that the cosmic-ray 
acceleration efficiency is ~ 0.2 for that Mach number. Since they defined cosmic-ray as 
particles with energies larger than those of thermal particles (or the injection energy for 
acceleration), the fraction of protons having energies of > rripC^ must be smaller than 0.2. 
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Although it is not certain whether equation (jl]) can be extrap olated down to th e injection 
energy, the adopted value of e = 0.05 is consistent with that of iRyu et al.l (120031 ) because it 
is smaller than 0.2. 

Fig. [1] shows the spectrum of nonthermal emission from accelerated particles at t = 
2x10^ yr. The distance to the model cluster is 78.4 Mpc (the distance to the Perseus clustei0). 
Synchrotron emission from primary electrons is dominant upto ~ 100 keV. The maximum 
energies for protons and electrons are -Emax.p 
respectively. The radius of the shock at this time is Rg 
Vo = 2650 km s~^, and the Mach number is 2.7. 



6.8 X 10^^ eV and E^axe = 2.1 x lO^^ ^y^ 



97 kpc, the shock velocity is 



Fig. [2] shows the spectrum at t = 4 x 10^ yr. At this time, 



max,p 



4.1 X 10^^ eV, 

-Smax,e = 2.1 X 10^^ cV, and Rg = 143 kpc, which is close to the size of the minihalo in 
the Perseus cluster. The shock velocity is Vg = 1950 km s~^ and the Mach number is 2.0. 
In Figs. [1] and [21 we plot the o bserved radio fluxes of the minihalo in the Perseus cluster 



( ISij bring 



1993 



Gitti et al.ll2002l ). As can be seen, the predicted radio synchrotron emission 



(long-dashed line) is too bright to be consistent with the observations. If we take smaller Eq, 
the radio luminosity becomes smaller. However, the size of the minihalo is too small to be 
consistent with the observed one. Moreover, if we consider a much smaller electron-proton 
ratio (e.g. re_p ~ 10~^), the radio spectral index is inconsistent with the observations. 

One possibility is that the age is much larger than 4 x 10'' yr and is t > 10^ yr. At 
that time, the shock is not prominent because it is almost in pressure equilibrium with the 
surrounding ICM. In fact, for the Perseus cluster, a shock of Rg ~ 100-200 kpc has not been 
reported. As the shock expands, its Mach number decreases from 2.7 (at t = 2 x 10'^ yr) 
to 2.0 (at t = 4 X 10"^ yr). Accordingly, the compression ratio (r^) decre ases, which affects 
the energy spectrum of particles and the emission from them (eq. [1]). iRyu et al.l (120031 ) 
indicated that particles are no longer accelerated if the Mach number is < 2. Thus, for 
t > 4 X 10'' yr, particle acceleration at the shock is not effective. 

At t ~ 10® yr, the emission from primary electrons may have died out, because electrons 
with a Lorentz factor of 7 > 10^, which are responsible for the radio emission, lose their 
energy t hrough synchr otron emission and inverse Compton emission on a time-scale of < 4 x 
10^ (e.g. ISarazinlll999l ). On the other hand, the lifetime of protons is much larger than 10® yr 
(eq. [7]), and the diffusion time of protons having the m aximum energy of ~ 10^^ eV from 
the c entral region of the cluster (~ 200 kpc) is 2 x 10® yr (jVolk. Aharonian. fc Breitschwerdt 



19961 ). Since particles are no longer accelerated at t > 4 x 10^ yr, the overall spectrum 



^The redshift of the Perseus cluster is 0.0183. We assumed that the cosmological parameters are fto — 0.3, 
Ao = 0.7, and Hq = 70 km s^^ Mpc^^ 
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originated from proton-proton collisions stays almost intact for 4 x 10^ ^ ^ ^ 2 x 10^ yr. 
In the radio band, only synchrotron emission from secondary electrons (dot-dashed line in 
Fig. [2]) will be observed at t ~ 10® yr. In this case, the predicted spectrum (thick-solid line 
in Fig. [2]) well fits the radio observations. Since the secondary radio emission is produced 
by protons having energies of ~ 100 GeV and since the diffusion time of these protons is 
> Gyr, the radio emission could persist that time. We emphasize that the assumption on the 
electron-proton ratio (^e-p) is not required to estimate the emission originated from proton 
collisions. 



4. Discussion 



Although our model is basically an one-zone model and cannot quantitatively predict 
the spatial change of the spectrum, we can qualitatively predict that. Compared with the 
one at t = 2 X 10^ yr (Fig.[T]), the spectrum of synchrotron emission from secondary electrons 
at t = 4 X 10'' yr is softer in the radio band (~ 0.3-1 GHz; Fig. [2]). The spectral index in the 
band of 327-609 MHz changes from 1.64 (t = 2 x 10^ yr, Rs = 97 kpc) to 1.88 (t = 4 x 10^ yr, 
Rs = 142 kpc). Since some of the protons accelerated at an earlier time should remain in 
the inner region of a cluster, the spectr um should be l ess steep in the inner region. This 
tendency is consistent with observations (jSijbringl Il993l : iGitti et al.l |2002| ) . 



As we mentioned above, the radio emission from secondary electrons could persist for 
a long time (> Gyr). Our mode l will be tested for the frequency (or the event rate) of 
large outbursts. iGitti et al.l (120071 ) indicated that large outbursts are likely occurring ~ 10% 
of the time in a significant proportion of all cooling core clusters. Thus, our model would 
suggest that a large fraction of clusters should have minihalos. This is inconsistent with the 
rareness of minihalos . However, the outbursts observed so far are of energies of < 10^^ erg 



( iRafferty et al.ll2006l ). which is smaller than our finding (1.8 x 10 erg). Thus, the rareness 



may indicate that strong AGN outbursts with energies of > 10^^ erg are rare phenomena 
or minor cluster mergers often perturb cluster cores. Another possibility is that particle 
acceleration at low-Mach number shocks occurs only in some specific IGM environments 
depending on the density of the surrounding matter, magnetic field configurations, and so 
on. In the future, statistical studies about AGN outbursts of ~ 10^^ erg are highly desired. 
The morphology of the radio surface brightness would also be important to check the validity 
of the model. If the particle acceleration is triggered by the expanding shock, one would 



expec t a torus-like shape instead of the spherical shape observed for minihalos (IGitti et al. 



2OO2I ). However, for clusters observed so far, the central region behind the shock is not 
empty; the IGM is still filling (e.g. Fig. 3 of iMcNamara et al.ll2005l ). Thermal protons there 
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may work as target protons for the proton-proton interaction and thus the radio emission 
may not be a torus-hke shape. A spatially resolved model must be constructed to address 
this issue. 

In Figs. [1] and [2], we also plot the obser vational upper limits of gamma-ray emission 



from the Perseus cluster (jPerkins et al.ll2006l ). At t = 2 x 10'' yr, the gamma-ray emission 



is brighter than the observations. At t ~ 10^ yr, there is no longer emission from primary 
electrons, and only gamma-ray emission of proton origin (thick-solid line in Fig. [2]) will 
be observed in the gamma-ray band. The predicted gamma-ray flux at ~ 10^ eV is 
~ 1 X 10~^^ erg cm~^ s~^, which could be detected with GL^S't[^ with a sensitivity of 
~ 3 X 10~^^ergcm~^s~^. The gamma-ray emission would persist for ~ tpp. If the gamma-ray 
is detected, it directly indicates that protons as well as electrons are accelerated in clusters. 
Moreover, the luminosity reflects the total energy of the protons. 

On the other hand, it would be difficult to detect the emission with imaging atmospheric 
Cherenkov tetecope. For example, H.E.S.sE has a sensitivity of ~ 1 x 10 erg cm ^ s ^ 
at ~ 10^^ eV). The predicted flux is smaller than the detection limit (Fig. [2]). 
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Fig. 1. — The spectrum of a shock at t = 2 x 10^ yr. Emissions from primary electrons 
are synchrotron (long-dashed), bremsstrahlung (short-dashed) and inverse Compton (thin 
dotted). Emissions related to protons are 7r°-decay gamma-ray (thin-solid), synchrotron 
(dot-dashed), bremsstrahlung (short-and-long dashed), and inverse Compton (thick-dotted) 
emissions from secondary electrons. Th e thick-solid li n e shows the tota l nonthermal flux. 
Radio observations are shown bv dots (ISiibrind Il993l : iGitti et al.l |2002| ). and gamma-ray 



upper limits are shown by arrows ( jPerkins et al.ll2006l ). 
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Fig. 2. — Same as Fig. [T], but for t = 4 x 10^ yr. In contrast with Fig. [H the thick-sohd hne 
shows the flux from protons. 



